Limits on the time variation of the electromagnetic fine-structure constant in the low 
energy limit from absorption lines in the spectra of distant quasars 
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Most of the successful physical theories rely on the constancy of few fundamental quantities (such 
as the speed of light, c, the fine-structure constant, a, the proton to electron mass ratio, fi, etc), and 
constraining the possible time variations of these fundamental quantities is an important step toward 
a complete physical theory. Time variation of a can be accurately probed using absorption lines seen 
in the spectra of distant quasars. Here, we present the results of a detailed many-multiplet analysis 
performed on a new sample of Mg ll systems observed in high quality quasar spectra obtained using 
the Very Large Telescope. The weighted mean value of the variation in a derived from our analysis 
over the redshift range 0.4 < z < 2.3 is Aa/a = (-0.06 ± 0.06) x 10"^. The median redshift of 
our sample (z~1.55) corresponds to a look-back time of 9.7 Gyr in the most favored cosmological 
model today. This gives a 3a limit, -2.5 x lO"^** yr"^ < (Aa/aAt) < -1-1.2 x 10"^" yr"^, for 
the time variation of a, that forms the strongest constraint obtained based on high redshift quasar 
absorption line systems. 

PACS numbers: 98.80.Es,06.20. Jr,98.62.Ra 



Contemporary theories of fundamental interactions al- 
low the fundamental constants to vary as a function of 
space and time 1]. Therefore, direct measurements of 
possible time variations of various fundamental quan- 
tities are of utmost importance for a complete under- 
standing of fundamental physics. In this letter we in- 
vestigate the time variation of the fine-structure con- 
stant (a= e^/hc = 1/137.03599976(50) 0, where e is the 
charge of the electron and h the reduced Planck constant 
measured in the laboratories on Earth). The time evolu- 
tion of a can be probed in the framework of standard Big- 
Bang models using measurements performed at different 
redshifts {z). The Oklo phenomenon |3j], a natural fission 
reactor that operated 2 Gyrs ago, or at z ~ 0.16, gives the 
strongest constraint on the variation of a over cosmolog- 
ical time-scales, [Aa/aAt] = (-0.2 ± 0.8) x 10-^^ yr^^. 
At higher redshifts, a possible time dependence will be 
registered in the form of small shifts in the absorption 
line spectra seen toward distant quasars Q as the en- 
ergy of the atomic transitions depend on a. Initial at- 
tempts to measure the variation of a were based on the 
absorption lines of alkali-doublets(5j(AD method). The 
best constraint obtained using this method is Aa/a = 
(-0.5 ± 1.3) x 10~^[i|. Other methods such as the one 
using O III emission lines 0, though more robust, are not 
sensitive enough to detect small variations in Aa/a . 
Investigations based on molecular lines detected in 
two systems give Aa/a = (—0.10 ± 0.22) x 10~^ and 
(-0.08 ±0.27) X 10-^ at Zabs= 0.2467 and 0.6847 respec- 
tively. Such studies at high-z are elusive due to lack of 
molecular systems. The generalization of alkali-doublet 
method, called many-multiplet (MM) method, gives an 
order of magnitude improvement in the measurement of 
Aa/a compared to AD method 9] by using not only dou- 



blets but several multiplets from different species. The 
sensitivity of different line transitions from different mul- 
tiplets to variations in a were computed using many- 
body calculations taking into account dominant relativis- 
tic effects [loj. 

In simple terms, MM method exploits the fact that 
the energy of different line transitions vary differently for 
a given change in a. For example rest wavelengths of 
Mg iiAA2797,2803 and Mg iA2852 transitions are fairly 
insensitive to small changes in a thereby providing good 
anchor for measuring the systemic redshift. Whereas the 
rest wavelengths of Fe ii multiplets are quite sensitive to 
small variations in a. Thus measuring consistent relative 
shifts between an anchor and different Fe ii lines can in 
principle lead to an accurate measurement of Aa/a. The 
accuracy at which the variation can be measured depends 
very much on how well absorption profiles can be mod- 
eled. For this usually Voigt profiles that are convolved 
with the instrumental profile and characterized by col- 
umn density {N), velocity dispersion (6) and redshift in 
addition to the rest-wavelength of the species are used. In 
a real spectrum small relative shifts can be introduced by 
various systematic effects such as inhomogeneities in the 
absorbing region, poor wavelength calibration, isotopic 
abundances, and atmospheric dispersion effects, etc[13|. 
However most of the random systematic effects can be 
canceled by using a large number of measurements. MM 
method applied to large samples of quasar absorption 
lines resulted in the claim for smaller value of a in the 
past^ij, Aa/a = (-0.574±0.102) x 10"^ for 0.2< z <3.7. 
The main motivation of the present work is to perform 
similar analysis using a completely different line fitting 
code and an independent, uniform, better quality, and 
well defined data sample. 
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The data used in this study were obtained with the 
Ultra-violet and Visible Echelle Spectrograph (UVES) 
mounted on the ESO Kueyen 8.2 m telescope at the 
Paranal observatory for the ESO-VLT Large Programme 
"QSO absorption lines". This programme has been de- 
vised to gather a homogeneous sample of echelle spectra 
of 18 quasars, with the uniform spectral coverage, resolu- 
tion and signal-to-noise ratio (S/N) suitable for studying 
the intergalactic medium in the redshift range 1.7—4.5. 
Spectra of quasars were obtained in service mode obser- 
vations spread over 4 periods (2 years) covering 30 nights 
under good seeing conditions (< 0.8 arcsec) using a slit 
1 arcsec wide. The data were reduced using the UVES 
pipeline, a standard set of procedures implemented in a 
dedicated context of MIDAS, the ESO data reduction 
package. However independent consistency checks were 
made to ensure good sky-background subtraction (better 
than 1% accuracy), accurate wavelength calibration (rms 
in AA/A better than ~ 4xl0~^ over the full wavelength 
range of interest, 3100-5400 and 5450-9000 A) and op- 
timum signal to noize (S/N ~ 50—80 per pixel). We fit 
the continuum using lower order polynomial that gives 
~ 1. Air- vacuum conversions and heliocentric correc- 
tions were done using standard conversion equations 0|. 
The slit was always oriented along the parallactic angle 
and calibration exposures were taken before or after the 
scientific exposures. The full descriptions of the data re- 
duction and calibration procedures, sample selection and 
all the steps involved in the analysis are presented in 
Chand et al., (2004) 12]. 

As a first step we have performed a detailed analysis on 
simulated data to have a clearer understanding of various 
possible systematics that can affect the analysis of real 
data. The results of this exercise are used to validate our 
procedure and define selection criteria that will minimise 
systematics in our analysis. Absorption spectra of Mg ii 
and Fe ii were simulated for given N, and b. The spec- 
tral resolution and S/N ratio were similar to our UVES 
data. The standard fitting function and coefficients^ 
are used to incorporate the effects of a variation. We con- 
sidered two cases: a simple single component system and 
a highly blended two-component system. In the highly 
blended case we restricted the separation between the 
two components to be always smaller than the velocity 
dispersion of one of the components. We then fitted the 
absorption lines in order to recover Aa/a introduced in 
the input spectrum. This exercise is to determine the 
precision that can be reached using our method and fit- 
ting procedure. For this we used the Voigt profile fitting 
method and standard statistics to fit the absorption 
profiles consistently and to determine the best fit value 
for Aa/a . The relationships between the input and re- 
covered values of Aa/a are shown in Fig.^ The method 
works very well in the case of simple single component 
systems where one expects minimum uncertainties due to 
systematics. The deviation of the recovered value with 
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FIG. 1: Absorption spectra of Mg ll and Fe ll were simu- 
lated for a given set of A'^, 6, and spectral resolution similar 
to our data, introducing spectral shifts corresponding to a 
given value of Aa/a . Top panels show the relationship be- 
tween the input and derived value of Aa/a in the case of a 
single clean component (left-hand side) and a blend of two 
components (right-hand side). Each realization is performed 
using random values oi N, b and noise keeping the signal-to- 
noise ratio, wavelength sampling and resolution as in a typical 
UVES spectrum. Dots are the values from individual realiza- 
tions and the points with the error bars are the weighted mean 
obtained from 30 realizations. A typical absorption profile is 
also shown in these panels. The lower panels give the distri- 
bution of the recovered Aa/a around the true one. Single 
(left) and blended (right) cases are considered respectively. 
Best fitted Gaussian distributions are over-plotted. 



respect to the true one distributes like a Gaussian with 
a = 0.23 X 10~^. This shows that Aa/a can be con- 
strained with an accuracy of 0.07xl0~^ even when we 
use 10 such systems. As expected the scatter is more in 
the case of blended systems {a = 0.38 x 10~^). Simi- 
lar analysis on two components with separations larger 
than individual b values gives a = 0.261 x 10^^. In- 
deed, in this case, the separation of the components is 
large enough so that each component is clearly recog- 
nized by the fitting procedure. In such a situation the 
corresponding errors are very much the same as for indi- 
vidual components. The conclusion is that it is better to 
select single component or not heavily blended systems, 
that is the components are separated by more than their 
6 values, to obtain the most reliable results. We notice 
from our simulations that systems with weak lines have 
large uncertainties in the measured Aa/a and should be 
avoided from the analysis. Finally, the data used in this 
analysis have a median S/N ratio of 70 per pixel, a factor 
of two better than that used in earlier studies. Based on 
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FIG. 2: Voigt profile fits to 4 randomly chosen systems (out 
of 23) in our sample are shown in the first 4 columns from the 
left. The quasar name and absorption redshifts are given on 
top of the panels. The points with the error bars are the ob- 
served data and the continuous curve is the fit obtained using 
multicomponent Voigt profile decomposition. The locations 
of different components are marked with vertical dotted lines. 
The plots in the right most column demonstrate how Aa/a is 
extracted from these systems. We plot as the function of 
Aof/a . The minimum in this curve (marked with a dot) gives 
the best- fit value of Aa/a and the error in this measurement 
(error-bar around the dot) is based on the standard statistical 
method of computing errors from A^^ = 1. 



the simulations we expect this enhancement in the S/N 
ratio to provide a roughly a factor two improvement in 
the Aa/a measurements. 

In summary, based on the results from the simulations, 
we apply the following selection criteria to derive reliable 
Aa/a : (i) we consider only species with similar ion- 
ization potentials (Mg ii, Fe il, Si ii and Al ll) as they 
are most likely to originate from similar regions in the 
cloud; (ii) we avoid absorption lines that are contami- 
nated by atmospheric lines; (iii) we consider only systems 
that have A^(Fe ii)> 2 x 10^^ cm~^ which ensures that all 
the standard Fe II multiplets are detected at more than 
5(7 level; (iv) we demand that at least one of the anchor 
lines is not saturated so that the redshift measurement is 
robust; (v) we also avoided sub-DLAs (i.e iV(H l)> 10^^ 
cm~^) as these systems may have ionization and chem- 
ical inhomogeneities; (vi) we do not consider strongly 
saturated systems with large velocity spread (complex 
blends) ; however in such systems whenever we find a well 
detached satellite components we include these compo- 



nents in the analysis; (vii) finally, based on the com- 
ponent structure resulting from the Voigt profile fits of 
systems that are not complex blends, we retain only sys- 
tems for which the majority of components are separated 
from its neighboring components by more than the b pa- 
rameters. Application of the above conditions resulted in 
23 systems for performing Aa/a measurements(6 single 
component systems, 6 well separated doubles, 6 systems 
with 3 components with at least one well detached from 
the rest, and 5 systems with more than 3 components) . 

In our analysis we use the fact that relative shifts be- 
tween lines from the same species (say Fe ii multiplets) 
are insensitive to Aa/a . We first fit all the Fe ii lines 
(except for Fe 11AI6O8) simultaneously using laboratory 
wavelengths. This allows us to find out about (i) bad 
pixels, (ii) unknown contaminations and (iii) the veloc- 
ity component structure in Fe 11. Defective absorption 
lines detected are removed from the Aa/a analysis. Sim- 
ilar exercise was carried out for Mg 11 doublets and other 
anchors. Based on these preliminary fits a first set of pa- 
rameters is generated to start the Voigt profile fitting pro- 
cedure that includes Aa/a variations. We illustrate the 
method in Fig. |21 using 4 randomly chosen systems from 
our sample. We vary Aa/a ranging from —5.0 x 10^^ to 
5.0 X 10~^ in step of 0.1 x 10^^, and each time fit all the 
lines together, minima obtained for each of these fits 
are plotted as a function of Aa/a (right most panel in 
Fig.|21). The value of Aa/a at which this is minimum 
is accepted as the measure of best possible Aa/a value. 
Following standard statistical procedure we assign Icr er- 
ror bar to the best fitted value of Aa/a by computing 
the change in Aa / a implying A^^ = ~ Xmin = 1 • We 
always consider two different models: (i) one in which h 
parameters for all the species are the same and (ii) one in 
which h parameters for different species are different. In 
all systems we notice that the derived Aa/a using both 
cases are consistent with one another within Icr uncer- 
tainty. We use the result that has a lower reduced for 
our final analysis. 

Results obtained for the 23 systems in our sample and 
that from the literature are summarized in Fig. |2| The 
shaded region passing through most of the error bars 
is the weighted mean (with 1/error^ weights) and its 
3(7 error from our sample. It is clear that most of our 
measurements are consistent with zero within the un- 
certainties. The simple mean, weighted mean and stan- 
dard deviation around the mean obtained for our sample 
are (-0.02 ± 0.10) x 10-^ (-0.06 ± 0.06) x 10"^ and 
0.41 xlO~^ respectively. The corresponding values are 
(0.01±0.15) X 10-^ (-0.08±0.07) x 10"^ and 0.26x10"^ 
respectively when we restrict our analysis to the single 
and well detached two component systems (12 in total). 
These systems have maximum infiuence in the weighted 
mean of the whole sample. Interestingly, error on the 
weighted mean and standard deviation around the mean 
from these systems are very close to the predictions from 
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FIG. 3: The measured values of Aa/a from our sample (filled 
circles) are plotted against the absorption redshifts of Mg ll 
systems. Each point is the best fitted value obtained for in- 
dividual systems using minimization as demonstrated in 
Fig. |5| The open circle and stars are the measurement from 
Oklo phenomenon^ and from molecular linesQ respectively. 
The weighted mean and la range measured by Murphy et 
al.(2003) are shown with the horizontal long dashed lines. 
Clearly most of our measurements are inconsistent with this 
range. The shadow region marks the weighted mean and its 
3(T error obtained from our study. 



our simulations. In all cases weighted mean values are 
consistent with all the data points with a reduced ~ 1. 

Thus our study gives a more stringent 3cr constraint 
of -0.24 < Aa/a{ in 10"^) < +0.12 over the redshift 
range of0.4< z <2.3. The median redshift of the whole 
sample is 1.55 which corresponds to a look-back time 
of 9.7 Gyr in the most favored cosmological model to- 
day (f7(total) = 1, f^A = 0.7, f7(matter) = 0.3, and Hq 
— 68 km s^^Mpc^^). This gives a 3cr constraint on 
the time variation of Aa/a to be —2.5 x 10~^^ yr^^ < 
(Aa/aAt) < 1.2 x lO"^*" yr"^ As can be seen from Fig. El 
our results are consistent with the constraints from Oklo 
phenomenon and from the molecular lines. However, our 
study does not support the claims by previous authors 
of a statistically significant change in Aa/a with cosmic 
time at z > 0.5 using MM method. In any case our mea- 
surement still does allow smaller variations in excess of 
what is found based on the Oklo phenomenon. Future 
very high resolution (R'^100,000) spectroscopic studies 



are needed to probe the variations in a with much better 
accuracy. 
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